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Abstract 

The critical current density shown by a superconductor at the extreme type-II hmit is predicted 
to follow a power law with external magnetic field B if the vortex lattice is weakly pinned 

by material line defects. It acquires an additional inverse dependence with thickness along the 
line direction once pinning of the interstitial vortex lines by material point defects is included. 
Moderate quantitative agreement with the critical current density shown by second-generation 
wires of high-temperature superconductors in kG magnetic fields is achieved at liquid-nitrogen 
temperature. 
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Thin films of tlie liigh-temperature superconductor YBa2Cu307_5 (YBCO) grown by 
pulsed laser deposition (PLD) on a substrate that serves to align the crystalline axes now 
routinely achieve critical currents that reach a substantial fraction of the maximum depairing 
. Such films also typically show a peak in the critical current as a function 
of the orientation of external magnetic field near the crystallographic c axis, perpendicular 
to the film (copper-oxygen) planejsl. It is well known that the Abrikosov vortex lattice 
that exists in external magnetic field shows no critical current if it is not pinned down by 
material defectsj^ The c-axis peak in the critical current shown by films of PLD- YBCO 
therefore indicates that some fraction of the material defects that pin the vortex lattice must 
be correlated along the c axis as well. Careful studies of the microstructure in such YBCO 
films find that the correlated pinning centers in question are in fact lines of dislocations 
th. appea. .a.u.a..v ,n *e PLD g™w.h .ro^eM- The pi,.i.g ce.te. no.abl. a„a„ge 
themselves in a manner that resembles a snapshot of a two-dimensional (2D) liquid, as 
opposed to a snapshot of a 2D gas. 

Below, we explore the theoretical consequences of the proposal that the vortex lattice 
induced by an external magnetic field oriented near the c axis in films of PLD- YBCO is in a 
thermodynamic Bose glass state characterized by an infinite tilt modulusjfl]. We specifically 
focus on the high-field regime at the extreme type-II limit, in which case only a fraction of the 
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vortex lines are localized at the dislocations that thread the film along the c axis 
and in which case the pinning of the vortex lattice is collective Q]. The theory predicts 
a critical current density along the film that approximately follows a power law as 

a function of magnetic field B in the collective pinning regime. It successfully accounts 
for the critical current density obtained in c-axis field from films of PLD- YBCO that are 
microns thick, at liquid nitrogen temperature. It fails, however, for much thinner films at 
lower temperature. We demonstrate that this failure can be corrected by including the effect 
of point pins on the interstitial vortex lines that lie in between the correlated pins. They 
contribute an inverse dependence on film thickness to the the critical current density in 
magnetic field oriented near the c axis. 

Consider again the film of YBCO just described in large enough external magnetic field 
applied along the c axis so that magnetic fluxlines overlap considerably. The vortex lines 
experience long range repulsive interactions in such case, which favors their arrangement 
into a triangular lattice. The natural material line defects that also run parallel to the c 
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axis act to pin down some fraction of the vortex lines, on the other hand. This frustrates 
the formation of a vortex lattice. The resulting thermodynamic state is a Bose glas s Ql that 
displays a divergent tilt modulus because of pinning by the material line defect s[10|jll| 
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The vortex lattice described above is effectively 2D due to the strong orientation along 
the c axis. Theoretical and numerical calculations then indicate that a net concentration of 
straight lines of unbound dislocations will be quenched into the triangular vortex lattice at 
zero temperature I 16 1. Monte Carlo simulations of the frustrated XY model further 
show that the dislocations in the vortex lattice appear either unbound or bound-up into 
neutral pairs in the extreme type-II limit Q]. In particular, the lines of dislocations do not 
arrange themselves into low-angle grain boundaries (cf. refs 17| and la])- 



We shall now compute the critical current density that is expected along the film direction 
for the hexatic Bose glass described above|ll| 3]. The vortex lattice is pinned collectively B 
when the number of vortex lines localized at material line defects is small compared to the 
total number of vortex lines. This shall be assumed throughout. The positional coherence 
of the vortex lattice along the c axis implies 2D collective pinning in particular BQ], with 
an infinite Larkin length in that direction 12]: Lc — oo. Next, we shall assume that the 
critical current density is limited by plastic creep of the vortex lattice due to slip of the 
quenched-in lines of dislocations along their respective glide planes The Larkin length 
in the direction transverse to the magnetic field. Re-, is then obtained by minimizing the sum 
of the elastic energy cost and the gain in pinning energy due to the translation of a Larkin 
domain by an elementary Burgers vector of the triangular vortex lattice, h = a/\. This yields 
the estimate lallldl 



CoUpifp/ceeby 



(1) 



for the density of Larkin domains, where Up denotes the density of pinned vortex lines, 
where fp denotes the maximum pinning force per unit length along a material line defect, 
and where Cqq denotes the elastic shear modulus. The prefactor Cq above is of order unity. 
On the other hand, because the Larkin correlation volume corresponds to the largest bundle 
of vortex lines that exhibits a purely elastic response ^ , it is natural to identify the average 
separation between lines of dislocations that are quenched into the vortex lattice with the 
transverse Larkin length, Rc. A similar minimization of the sum of the elastic energy and 
the pinning energy then yields the same form (P) for the density of Larkin domains in such 
casejlj], but with a prefactor Cq = 7r/ln(i?c/<^df)^- H^re a^f is of order the core diameter 



of a dislocation in the vortex lattice. Finally, the statistical nature of 2D collective pinning 
requires many pinned vortex lines per Larkin domain : Up > R^"^ . Comparison with 
relation (Q) then yields the threshold in magnetic field 

5cp = C2(v^/2)(4/p/£o)'$o, (2) 

above which 2D collective pinning holds. The estimate Cge = ($o/8^'^l)^'^vx for the shear 
modulus of the vortex lattice [2^ has been used here, where rivx denotes the density of vortex 
lines, where \l denotes the London penetration depth, and where $0 denotes the flux quan- 
tum. Above, = {^al '^t^^lY is the maximum tension of a fiuxline in the superconductor. 
The Lorentz force balances the pinning force in the critical state following jcB/c = 
■ ripY^'^fp when 2D collective pinning holdspl]. Here jc denotes the critical current 
density along the film, which is perpendicular to the magnetic field B aligned parallel to 
the material line defects. Substitution of the result for the density of Larkin domains 
reduces this balance to 

j,B/c = Coiceeby'ripfl (3) 

The critical current density can then be determined from the above condition for mechanical 
equilibrium once the density of pinned vortex lines, Up, is known. Non-interacting vortex 
lines, Cee/Zp — > 0, yields the upper bound for Up. The profile of Up versus magnetic field in 
this case is clearly just the upper pair of dashed lines that join at the density of material line 
defects, ra^, shown in Fig. [T] Infinitely weak correlated pins, Cge/Zp 00, yields the lower 
bound for Up, on the other hand. It corresponds to the limit of an infinitely rigid vortex 
lattice. The profile of np versus magnetic field in that case follows the lower dashed line 
shown in Fig. ^ where Vp denotes the effective pinning radius of the material line defects. 
It is obtained by observing that Up/n^ coincides with the probability that a given pinning 
center trap a vortex line, p = vrrp ■ n^^. (See Fig. |21and ref. ||2l|.) This coincidence is true 
for any density of pinned vortex lines. Hp, as long as the material line defects do not crowd 
together: nr^ ■ <^ 1. Such is the case with the lines of dislocations that thread films of 
PLD-YBCO, which display an effective hard-core repulsionj^] ^. 

The fact that the vortex lattice is well ordered within a Larkin domain (see Fig. El) indi- 
cates that the direct scaling of Hp with B obtained above in the limit of an infinitely rigid vor- 
tex lattice persists in the regime of weak correlated pinning more generally. This in fact can 
be demonstrated by first observing that the relative statistical error in the number of pinned 



vortex lines at /p — > obeys the law of large numbers: AN^/Np = /(vrrp ■ n^x)/ \/n^Rci 
where Np denotes the number of pinned vortex lines inside of a Larkin domain, and 
where f{p) = J (1 —p)/p sets the fluctuation scale specifically for the binomial probability 



distribution |2^. Second, adapting the statistical principle of Larkin-Ovchinnikov[6|] to the 
present case yields the new rule that the break-up of the pristine vortex lattice into Larkin 
domains of dimensions Rc x Rc serves to increase the number of pinned vortex lines with 
respect to the most probable number, Np = p ■ n^R^^. at /p ^ 0, by a number of order the 



statistical error, AA^pi 



Here Cq is a numerical constant of order unity. The density of pinned vortex lines therefore 



scales nearly directly with magnetic field at small relative errors, ANp/Np <^ 1. This 
requires collective pinning, Np ^ 1. Substitution into Eq. yields a density of Larkin 
domains, R^^ = (-\/37r/2)CQ (4/prp/eo)^'T-(/., that depends only weakly on the magnetic field 
in such case. Finally, the function f{p) specific to the binomial probability distribution can 
be approximated by (2/v^) ln[{3^/e/A)/p] for probabilities p inside of the range [0.01,0.98] 
(cf. ref. |2l|). Substituting that approximation into Eq. (jlj) then yields the scaling result 

rip/n^ = {B/B^^Y^'^ (5) 
for the density of pinned vortex lines as a function of magnetic field oriented parallel to 



the material line defects, where dp/2 = 1 — (ci/y ra^i?^) gives one-half the effective scaling 
dimension, and where -002 = C2$o/7i"^p gives the saturation field. Here, we define Ci = 
2co/v^ and C2 = (3v^/4)"['^^/'^p)^^l. The last constant is notably less than unity. Fig. [T] 
depicts the field dependence of the scaling ansatz Q schematically, where the correction to 
the scaling dimension is ignored. 

Table E] lists various limiting cases that obey the scaling ansatz for the density of 
pinned vortex lines as a function of magnetic field, including the lower limit at weak pinning 
just treated. Substituting (0) into the condition (jS)) for the critical state then yields the 
power law 

Je(S)/j,(5o) = (5o/i?)"- (6) 

for the critical current density as a function of magnetic field, with exponent 

«cp = (3 - cip)/2. (7) 



The reference critical current density jciBo) at the geometric mean Bq = (-Bcp-B<^2)^^^ of the 
range in magnetic field of the power law, [-Bcp, -B</,2], determines the latter by the formula 

V ' VBJT72' UBo)) ■ ^ ^ 

Here jo = 4c£:o/3-\/3$o^ and ifc2 = '^'o/Stt^^ are the the depairing current density and the 
upper critical field, respectively, each set by the coherence length ^ jTj]. Also, the constant 
factor above is defined by Ci = {16^/^/3^/^)Co. Study of Eqs. (0) and again shows that 
Up scales nearly directly with B in the collective pinning regime, B > Bcp. By the previous 
Eqs. (jSI)-©, we conclude that the critical current density then decays with magnetic field 
like 1//B (see Table P). 

Consider next the addition of a field of material point defects to the thin-film super- 
conductor. Either the Bose glass state described above remains intact 91, in which case 

n 

Lc remains infinite |12l|. or it will break up into smaller Larkin domains of finite thickness, 
Lc < 00. In either case, consider henceforth films of thickness r < Lc- Collective pinning 
therefore remains 2Dj^^^, and the interstitial vortex lines that remain free of the material 
line defects may be considered as rigid rods. The material point defects will pin the latter. 
Recall now the basic idea behind 2D collective pinning, which is that the critical pinning 
force per unit volume is given by the quotient of the variance of the net equilibrium force per 
unit length over a Larkin domain with the cross- sectional area R^. The magnitude-square 
of this quotient is equal to 



I E m\VRt = in,fp+n'X)/Rl (9) 

Above, / denotes the force per unit length due either to a vortex line pinned by a material 
line defect or to an interstitial vortex line pinned by material point defects. These have 
respective root-mean-square values of /p and /p. Also, and n'^ denote, respectively, the 
density of vortex lines pinned by material line defects and the density of interstitial vortex 
lines pinned by material point defects. Last, the overbar notation above denotes a bulk 
average achieved by rigid translations of a given Larkin domain (LD). The identity is 
then due to the statistical independence of all of the pinning forces: Z^ijgLD /(O ' fU) = 0) 
where the prime over the summation symbol specifies that i and j denote distinct pinned 
vortex lines within a given Larkin domain. This last property is closely related to the null 
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average value of the net pinning force exerted on a Larkin domain when the vortex lattice 
is in mechanical equilibrium. In particular, / = implies I^ld / = 0- By Eq. Q, the 
additional source of pinning due to the interstitial vortex lines can then be accounted for by 
making the replacements 

^p/p ^ ^p/p + n'pfp and jc Jc + j'c (10) 

in Eqs. and (jH)) for the density of Larkin domains in the vortex lattice and for the critical 
state. Here, j'^ denotes the contribution by the point pins to the net critical current density, 
Jc = jc + j'c- The first replacement displayed by Eq. (fTUj) indicates an effective density of 
vortex lines pinned by material line defects, + {f^/ /p)^p- Notice that the latter correctly 
equals when f'^ = 0, and Up + rip when f^ = /p. Matching it to the modified result for the 
density of Larkin domains ((H) then yields the old result for the threshold field beyond 
which 2D collective-pinning holds. The previous Eqs. (jHl) and (jHl) then remain unchanged 
as long as jc is understood to represent the contribution to the critical current density by 
the material line defects alone! 

The condition (jS)) for the critical state that results from pinning of the vortex lattice 
by material line defects demonstrates that the critical current density jc is a bulk quantity 
that is independent of the thickness of the thin-film superconductor along the axis of the 
correlated pins. This notably is not the case for the contribution by point pins to the critical 
current density, j'^. In particular, the forces due to point pins add up statistically along a 
rigid interstitial vortex line if the film is much thicker than the average separation between 
such pins along the field direction, t^. The effective pinning force per unit length experienced 
by an interstitial vortex line is then given byjl^ /p = fUi'^p'^Y^'^ at film thicknesses r > Tp, 
where /g denotes the maximum force exerted by a point pin. Also, all of the interstitial 
vortex lines are pinned by point defects if the film is thick enough: n'^ = n^^ — np at r > Tp. 
The relative contribution by point pins to the critical current density is then predicted 
to show an inverse dependence on film thickness, j'c/jc = tq/'^, that is set by the scale 
''o = {^vx/np — l){fQ / fpT^). The net critical current density, Jc = jc+j'c then follows a pure 
1 / \^ power law in the case that the density of vortex lines pinned by material line defects 
scales directly with magnetic field [dp = 2). More generally, we predict a linear dependence 
on film thickness for the net critical current per unit width following Kc = {tq + r)jc, where 
Kc = tJc by definition. 
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The summation of pinning forces is coherent, on the other hand, at magnetic fields below 
the threshold for 2D collective pinning: jcB/c = n^fp at B < B^p. Further, we have that 
Hp is approximately equal to n^^ at low magnetic fields compared to the accommodation 
scale[9] B^i = (4£:p/£:o)('^'o " ^(t>), where Ep denotes the depth of the correlated pinning 
potential per unit length. The above balance of forces then yields a critical current density 
that reaches a plateau 

Jc(0+)/jo = (3v^/4)(/pe/£o) (11) 



in the zero- field limit j23|]. Finally, the inequality B^i < B^p shall be assumed throughout. 
It places the bound -B^ < (4£:p/£:o)('^'o/''^p) on the matching field. 

Thin films of YBCO grown by PLD on a substrate also show an extended regime in 
magnetic field oriented parallel to the c axis where the critical current obeys a power law 
(see Fig. Eland refs. jll^PSl)- The films themselves are divided into columns of subgrains 
that run parallel to the crystallographic c axis. Etching of the film surface demonstrates 
that lines of dislocations also run parallel to the c axis in "trenches" that separate such 
growth islands^ There typically exists about one threading dislocation per growth 
island, each separated by a distance of about 110 nm Qlj^- The dimensionality of the 
columnar pins along a cross section then coincides with that of the subgrains, dp = 2. The 
critical current that is obtained from such films of YBCO in magnetic field B aligned along 
the c-a.,s typically obeys a l/VB law. i„ pa.t,c..a.QaaQ. Th,s is consistent with a 
direct dependence, Up oc B, of the density of pinned vortex lines on the magnetic field by 
the present theory (0) for 2D collective pinning. Such a dependence on magnetic field, in 
turn, is predicted approximately by Eq. (0)) in the collective pinning regime. (See Table EJ. 

Another point of comparison between theory and experiment is the actual magnitude of 
the critical current density. The maximum force per unit length exerted by a correlated pin, 
/p, can be obtained from jc in self field via Eq. (fTTjl |24| . Substituting it into the right-hand 
side of Eq. ^ then yields a prediction for the threshold field beyond which 2D collective 
pinning holds: 

Bcp = C',[UO+)/jo]'H,2. (12) 

Values of self-field Jc measured in various films of PLD- YBCO are listed in Table 

HTl along side of the predicted threshold field, B^p. The London penetration depth is set to 
its value at zero temperature, Al(0) = 150 nm. A value for the superconducting coherence 
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length of ^(0) = 1.5 nm then yields a depairing current density of jo(0) = 300 MA/cm^, 
as well as an upper-critical field of i/c2(0) = 146 T. Also, the prefactor Co on the right- 
hand side of Eq. (Q) for the density of Larkin domains in the vortex lattice is set to unity. 



The predicted threshold magnetic field fails to lie below the geometric mean y/BjBo of the 




power-law regime (see Fig. ^ only in the case of the relatively thin YBCO filmj2]j2i|^- This 
is likely due to the neglect of pinning by true grain boundaries 3] • 

The present theory for 2D collective pinning by material line defects also predicts the 
range of the power law in magnetic field (jH)) from knowledge of the critical current density 
at the logarithmic midpoint. (See Fig. IHl "x".) Table ITTl again lists recent measurements 
of the critical current density in films of PLD-YBCO at the logarithmic midpoint, along 
side of the predicted range. Only films that show an inverse power-law (JH)) in Jc versus B 
characterized by an exponent a in the vicinity of 1/2 are considered. Standard physical 
parameters for YBCO are again used, the prefactor Cq is again set to unity, and Eq is again 
set to its value at zero temperature. Additionally, the effect of point pins is neglected 
(jc — 0)- Comparison of the last two columns in Table ITTl indicates that the present theory 
for 2D collective pinning gives a fair account of the critical current density in PLD films of 
YBCO superconductor that are microns thick. We believe that the failure of the theory in 
the case of the much thinner film is due to the neglect of the contribution by point pins. 

In conclusion, 2D collective pinning of the vortex lattice by material line defects can 
account for the inverse-square-root power law obeyed by the critical current density in films 
of PLD-YBCO versus external magnetic fieldQ. We also predict on this basis that the 
critical current per unit width is a linear function of film thickness, with a positive slope 
equal to the bulk critical current density, and with a negative intercept on the thickness 
axis, after extrapolation, due to pinning of interstitial vortex lines by material point defects. 
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FIG. 3: Sketched is a typical profile of the critical current density measured in films of PLD-YBCO 
as a function of external magnetic field oriented near the c axis. The "x" marks the power-law 
regime in magnetic field, [Bi,B2]. 
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TABLE I: Listed are various examples of pinning by material line defects that obey the scaling 
ansatz Q for the density of pinned vortex lines as a function of magnetic field. 
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"Rough estimate taken from ref. 



TABLE II: Listed above are theoretical predictions, [i?cp, -6,^2]) versus experiment, [Bi,B2]. (See 
Fig. 121) The effect of point pinning is neglected {j'^ = 0). 



13 



